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Research on deformations of laser-welded 
joint of a steel sandwich structure model 


Dariusz Boronski 

Academy of Agriculture & Engineering, Bydgoszcz 
Janusz Kozak 

Gdańsk University of Technology 


ABSTRACT 


Problems of the behaviour of all-steel, laser welded sandwich panels in plate-to-web contact zone are 
presented. Results of laboratory tests on determination of deformation field around laser weld under ben- 
ding load — obtained by laser interferometry technique — are shown. Analysis of strain field distribution is 
included and some conclusions related to modeling of behaviour of such connection are drawn. 


Key words : laser welding, laboratory tests of structural elements, laser interferometry 


INTRODUCTION 


Developments in novel joining techniques including laser, 
hybrid (laser — MAG) and friction stir welding have brought 
about their broader and broader industrial applications, as well 
as they opened new possibilities for designers to create struc- 
tures not feasible so far. 

Higher accuracy of so made welds and their more uniform 
quality, increased speed of welding, reduced amount of heat 
introduced to weld zone, possible automation and use of ro- 
bots for welding processes — these are some advantages deci- 
ding upon fast application of such techniques in space, aircraft, 
car or shipbuilding industries [1, 2, 3, 4, 5, 6]. 

Their special position in the case of the application e.g. in 
shipbuilding [7] or aircraft industry is associated with, a.o., the 
possibility of building the sandwich structures which consist 
of load-carrying shell platings stiffened by cores of various 
forms, that leads, a.o., to a significant reduction of mass of 
designed objects. However application of the novel solutions 
usually demands some changes to be introduced to the tradi- 
tional methods of construction of such objects as aircrafts and 
ships. 

One of the example solutions based on the laser welding 
technique is a steel double-skin structure stiffened by internal 
webs joined with the shell platings by using laser welding ap- 
plied from outside of the object under construction. 

The idea of the replacement of the traditional single-skin 
ship hull structure with a new thin double-skin structure whose 
system of basic stiffening members is installed in the space 
between the shell platings, was born in the 1950s, and then it 
was studied by NASA [8, 9, 10]; however its practical applica- 
tion was attempted by US Navy at the end of the 1980s [11]. 
Presently Meyer Werft, Germany, applied the solution to the 
structural panels which appeared ten times stiffer and weighing 
less by 35% than their conventional equivalents [12]. 


Such panels have been produced and applied on an indu- 
strial scale e.g. to decks of Danube navigation passenger ships 
built in series by a German shipyard, as well as to bulkhead 
structures of another passenger ship [13]. 


a) b) 


12.5 mm 


200 x 90 x 9 mm 


Fig.1. Comparison of a classical deck structure (a) of a ro-ro ship 
and that built of Meyer 5 laser-welded panels (b) 


However in order to apply the novel structural elements on 
a large scale — especially to important elements of ship hull 
structure — it is necessary to have at one’s disposal information 
about their corrosion and fire resistance, and first of all about 
their technological and strength properties required to obtain 
approval of ship classification institutions. 

The problems have become the theme of the „Sandwich” 
research program being realized in the scope of 5» Frame Pro- 
gram financially supported by European Commission. Within 
the research program a team of the Faculty of Ocean Engine- 
ering and Ship Technology, Gdańsk University of Technology, 
carried out verification of strength properties of full-scale do- 
uble-skin panels internally stiffened by real stiffening member 
systems and different types of fillers [14]. Some new research 
problems which have resulted from the above mentioned work, 
are presently continued in the frame of the ASPIS project re- 
alized by the Faculty within EUREKA E!3074 program. 

The steel double-skin panels differ from their classical sin- 
gle-skin counterparts by their shear stiffness in the direction of 
the stiffeners (webs), much greater than that in the direction 
perpendicular to them, i.e. Dqgx>>Dqy. Simultaneously the stif- 
fness is much smaller than that of the relevant classical struc- 
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ture. It means that apart from the collapse mechanisms which 
may occur in one of the stiffened platings, similar to those po- 
ssible to occur in the traditional single-skin structures, other 
mechanisms resulting from mutual interaction of both platings 
may additionally occur, e.g. a kind of ,,rotation” of the ele- 
ments relative to each other around the laser weld as a,,hinge” 


(Fig.2). 


Beginning Process 
of contact of indentity 


mm ie ie | A 


Fig. 2. Behaviour of laser-welded joint under transverse load 


Bending of weld 


The phenomenon may be very important for the modelling 
of behaviour of a structure by means of the finite element me- 
thod (FEM). By neglecting it an excessively high stiffness cha- 
racteristics of a modeled structure against that real, may be 
erroneously calculated. Therefore from the point of view of mo- 
delling the behaviour of the entire panel under load it is impor- 
tant to know the macro-scale response of the joint. 

The laser-welded joint much differs from a conventional 
arc-welded joint. As a result of the laser welding a narrow weld 
surrounded by a small heat-affected zone but with entirely dif- 
ferent material properties against those of native material, is 
obtained. 

No universal guidelines (including fatigue test data) have 
been so far available for fatigue assessment of laser welds, as 
this is the case for the traditionally made welds [15]. The avai- 
lable results of high-cycle fatigue tests of unnotched specimens 
containing laser weld [16] indicate that their fatigue strength 
does not much differ from that of the native material, which 
can partially result from the better quality (lack of defects) of 
the weld itself. The result can be also explained on the basis of 
the local strain approach to fatigue strength in the elastic ran- 
ge. However in the case of the load causing exceedance of the 
material yield point in the joint, strain values in particular weld 
zones will be different. It should be also taken into account 
that plastic strains may occurr not only due to large external 
loads but also as a result of a complex geometry of the joint 
leading to strain/stress concentration. Obviously all the above 
mentioned effects may occur in the considered steel sandwich 
panels. 

In the reported work, being a part of the broader research 
program, the results of the first task of the research on local 
strain distribution in weld zone, due to the in-plane loading 
applied to the shell platings of the sandwich model structure, 
are presented. In the next phase of the investigations it is inten- 
ded to analyze strains developed in weld zone during fatigue 
tests. 


OBJECT AND METHOD 
OF THE RESEARCH 


The research program was aimed at revealing the local phe- 
nomena in the laser-welded joint under bending in order to 
collect information necessary for modelling the joint’s beha- 
viour. The object (specimen) used in the research represented 
a fragment of a steel sandwich panel. In Fig.3 its structural 
arrangement and basic dimensions are shown. In its measure- 
ment part the specimen contained three webs joined with both 
shell platings. The investigations were focused on the joint 
connecting the middle web with one of the shell platings. The 
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specimens were made of hull structural steel of known mecha- 
nical properties. 


a) 


80 mm 80 mm 


20mm 4mm 3mm 


is 
5 

e 
B 


Fig. 3. The structural arrangement (a) and the main dimensions (b) 
of the specimens 


The laser grating interferometry technique [17, 18] applied 
in the automated laser grating extensometer LES [19, 20], was 
used to investigate strains in the joint and to account for local 
character of the phenomena occurring in a small area of the 
joint. 

The technique is an optical full-field method based on the 
phenomenon of interference of two mutually coherent light 
beams projected on a measurement grating, which produce 
a pattern of interference lines (fringes) carrying information 
on relative displacements of the surface of a tested object. Sen- 
sitivity of the method depends on an applied laser wave-length 
and diffraction grating density. In the case of the LES it amo- 
unts to 0.42 um per one interference line. Further analysis of 
the lines makes it possible to obtain the measurement sensiti- 
vity of the order of 20 nm. 

The used research instrumentation consisted of the LES ex- 
tensometer, a fatigue testing machine and PC- class computer 
together with its accompanying equipment, is shown in Fig. 4. 


e PC with PCI boards: 
- NI PCI-GPIB 
- MVDelta (MVSigma) 
- Advantech PC 1750 


e fatique machine with 
digital control system 


e measurement head 


© software 
y © videorecorder 
and TV monitor 


e drivers and 
supply units for : 
- stepping motors 
- laser diodes 
- CCD camera 


Fig. 4. Configuration of the used research instrumentation 


Necessary calibration of the optical system can be remote- 
ly performed within the LES head. The pattern of interference 
lines (fringes) obtained during testing can be saved — through 
the electronic-optical vision system — in the computer memo- 
ry. The pattern is partially analyzed in the real-time mode, which 
makes it possible to measure strains in a selected segment of 
measurement field, and the full information is saved on ma- 
gnetic disk for further off-line analysis. Setting and operation 


of the entire research instrumentation can be controlled by 
means of a special software. 

The using of the laser interferometry method demands ap- 
plication of gratings placed onto the specimen’s surface. In the 
investigations in question the gratings of the density of 1200 
lines/mm obtained by dusting aluminium onto a matrix prepa- 
red on 2’ x 2’ glass plate, were used. During the preparation of 
the specimens for testing the 8 mm x 8 mm fragments of the 
grating were put on the cleaned surface of the joint by glueing 
the dusted aluminium layer, and then separating it from the 
matrix after setting the resin. 


TEST RESULTS AND THEIR ANALYSIS 


During the performed tests the distributions of relative dis- 
placements and strains in two mutually perpendicular direc- 
tions U and V for 16 selected levels of static load (complying 
with Fig.6 and 8), were determined for values of the loading 
being under control. The tested specimen fitted with the LES 
head is shown in Fig.5. The specimens were loaded in accor- 
dance with the loading scheme shown in Fig.6. 


ery 


Fig. 5. The sandwich structure model (specimen) under testing 


s] 


Loading levels 


Fig. 6. Loading scheme of the specimen 


At particular values of the loading force P the fringe pat- 
terns in two analyzed directions U and V were recorded. The 
measurements were carried out for the middle joint only. The 
strain measurement field is shown in Fig.7. 


a) b) 
P 


Investigated 
area 


Grid 


P 


Fig. 7. The specimen with the glued laser grating (a), 
and the measurement field location relative to the laser weld (b) 


The diagram of the load applied to a specimen in function 
of its elongation, together with the marked points of strain 
measurement, is presented in Fig.8. In the diagram the perma- 
nent elongation which remained in the specimen after unlo- 
ading, is additionally indicated. 


4.5 


4 i 


3.5 As 


N 
N wn U 
> 


Load [kN] 


0.2 0.4 0.6 0.8 1 
Elongation [mm] 
Fig.8. The diagram of the load applied 
to a specimen in function of its elongation 


The example fringe patterns recorded during the tests 


it HHN i i i i 


IAG i! 


U-direction 
V-direction 


P=2.5 kN 


A F f 
Fig.9. The example fringe patterns recorde: 
at P = 0 kN, and at P = 2.5 kN 
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The strain analysis was generally limited to the values of 
the force P = 2.5 KN as large strain values occurred in the con- 
tact zone of the web and shell plate, which exceeded the basic 
measuring range of the applied method (for the mode of auto- 
matic measuring and data recording), as well as due to the con- 
centration of fringes resulting from rotation of the specimen’s 
elements under loading. 

The obtained fringe patterns were analyzed with the use of 
numerical procedures making it possible to identify particular 
fringes and their phases within an analyzed field. 

In Fig.10 and 11 are presented the strain distributions de- 
termined in the directions V and U, respectively, for successi- 
ve values of the external force P. 


P=0kN P=0.23 kN 


P=0.55 kN 
oe 


The largest values of tensile strains (along V- axis) in the 
analyzed area occurred at the edge of the weld, and those 
of compressive strains — in the web-to-plate contact zone 
(Fig.12). 

Along with increasing the load the increasing strain values 
occurred both in the compression and tension zones, also 
a gradual expansion of the zones into the plate material was 
observed. 

The combined influence of the web bending and the web- 
-to-plate contact on the strain distribution in the minimum cross- 
-section of the weld, is shown in Fig. 13. 

From the analysis of the presented strain distribution it re- 
sults that the neutral axis of bending, located close to the middle 


P = 0.78 kN P = 0.98 kN 


y 


- 0.55 % 


Fig.10. Strain distributions in V-direction 


P=0.55 kN 


ere 


P = 0.78 kN 


J i sy? ae 
EJE: 


£] 


Fig.11. Strain distributions in U-direction 


From the analysis of the obtained strain distributions along 
V-axis it can be observed that the axial load applied to the 
specimen in accordance with the loading scheme (Fig.8), ge- 
nerated a strain gradient typical for bending, in the surroun- 
dings of the weld joining the web with the shell plate.Simulta- 
neously the edge of the web exerted contact pressure on the 
shell plate, that produced the compressive strain zones appe- 
aring both in the web and shell plate in the vicinity of the 
web’s edge. (Fig.12). 
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of the weld, did not change its location despite the web came in 
contact with the plate during loading the specimen. As a result, 
only a deviation from linear distribution in the compressive 
strain zone of the weld occurred. Simultaneuously the contact 
pressure of the web on the shell plate caused large plastic stra- 
ins in the contact zone, which exceeded the strains in the weld 
itself. Hence, such development of strains may indirectly indi- 
cate that the yield points of the weld, heat affected zone and 
native material could significantly differ from each other. 


Approximate 
position of weld 


Maximum tension 
strain zone 


Maximum compression 
strain zone 


Fig.12. Strain distribution in the weld zone at P = 1.54 kN 


0.3 
0.2 -r 
ON 
i, Weld Contact zone 
Ny, 
My, 
“ty 
0 “My 
Mays 
MIM a 
E€ Ny 
5-0.1 7) 
E D 
z t 
i—i ‘4, 
w yy 
g -0.2 4, 
3 4 
È LA 
NM -0.3 4 
“» 
Y 
ry 
0.4 th 
: -2 
4 
Y 
) 
-0.5 ha, i 
CAN 
0.6 
0.5 1 LS 2 
x [mm] 


Fig. 13. The strain distribution in the minimum 
cross-section of the weld, at P = 1.54 kN 

Due to lack of a firm connection between the structural 
elements in the contact zone the plastic strains in the native 
material, resulting from compression remain unchanged, how- 
ever the strains in the weld itself become smaller, that additio- 
nally increases the strain concentration effect in the heat affec- 
ted zone during unloading the specimen. 

The analysis of the strains in U-direction shows that they 
are symmetrically distributed with respect to the line run- 


a) 


— 


0.06 % 


y ON | _-— Laser weld 


Shell plating 


ning through the weld cross-section of the minimum area 
(Fig.14). 

After detail analysis of the typical strain distribution 
shown in Fig. 14.a, it was possible to suggest a probable cha- 
racter of the deformation of the joint containing laser weld, 
which is schematically shown in Fig.14.b. The web, while 
displacing relative to the shell plating, makes the laser weld 
formed of the material hardened due to welding process, ro- 
tating. The weld — while rotating — develops bending of the 
plate, manifested by two zones : of compressive and tensile 
strains, respectively, in the plate and web material on both 
sides of the weld. As a result, in the heat affected zone the 
strain gradients similar to those in the case of the strains in 
V-direction, appear. 


CONCLUSIONS 


> The investigation method applied in this work makes it 
possible to effectively determine strain distributions in the 
specimens of steel sandwich panels. The limitations asso- 
ciated with the exceedance of the basic measuring range of 
the method can be removed by using a sequential measure- 
ment procedure for the web and shell plate separately. 

> The obtained strain distributions revealed much lower ben- 
ding stiffness of the laser weld against the surrounding ma- 
terial, that made the entire joint more prone to deform. Due 
to a mutual ,,rotation” of the joined elements their edges 
(boundaries) come into mechanical contact in a short time. 

> The contact phenomena revealed during the tests and the 
resulting asymmetrical behaviour of the joint imply to be 
very cautious in using the finite element method for mo- 
delling structural behaviour of the objects containing such 
joints. Application of the rutine approach and typical finite 
elements for structural analysis of such objects may lead to 
erroneous, and even unsafe results. 

> In order to relate the observed effects associated with the 
developing of strain gradients to estimation of fatigue life 
of welded joints it would be necessary to perform a com- 
prehensive fatigue test program covering, a.o., determina- 
tion of weld-zone material properties under cyclic load and 
their changes during the testing. An analysis of the test re- 
sults would make it possible to verify the present approach 
to the traditional welded joints. Knowing the local state of 
strains in the weld one would be able to apply the strain or 
energy approach to calculation of fatigue life of structures. 


b) 
Web man 


Fig.14. Strain distribution in U-direction (a), and scheme of deformation of the joint (b) 
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ACRONIMS 


ASPIS - Application of Steel Panels into Ship Structure 
EMAS - Publisher 

FEM - Finite element method 

FPSOs - Floating Production, Storage and Offloading Units 
IW _ - International Institute of Welding 

IRCN - Institut de Recherches de la Construction Navale 
LES - Laser Grid Extensometer 

MAG - Metal - Active - Gas 

NASA - National Aeronautics & Space Agency 

RINA - Royal Intitution of Naval Architects 

SEM Society for Experimental Mechanics 

TWI - The Welding Institute 

VTT Valtion Teknillinen Tutkimuskeskus 
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Maritime Transport 2003 


On 25+28 November 2003 in Barcelona 
2° International Conference on : 


Maritime Transport & Maritime History 


organized by Department of Nautical Science 
and Engineering, Technical University of Catalonia, 
and Maritime Museum in Barcelona, had place. 


During the Conference 72 papers were presented. 


Most of the papers (14 and 9) concerned maritime trans- 
port and environmental protection, whereas 5 presenta- 
tions reminding maritime history which allows understan- 
ding ways of development and contemporary state of the 
art of the maritime knowledge and activities, deserved 
special attention. The presentation was enriched by visi- 
ting the Maritime Museum. 

Two sessions were carried out in Tarragona where the 
Conference participants had the occasion to visit the local 
Mediterranean port tightly connected with the important 
Spanish petrochemical complex and industrial centre. 

Prof. W. Galor of Maritime University of Szczecin (Po- 
land) took part in the activity of International Scientific 
Committee of the Conference; he also chaired one of the 
sessions on "Simulators and control systems" and presen- 
ted - during the session on "Training and simulations" - 
the paper on : "The navigational analysis of moderniza- 
tion of Leba port entrance". 

Moreover, during the session on "Ship development 
and hydrodynamics" dr. Jozef Kozak, Gdansk University 
of Technology (Poland), read the paper on : "Strength te- 


\ sts of steel sandwich panels". 


MARINE ENGINEERING 


Viscoelastic unsteady lubrication of radial slide 
journal bearing at impulsive motion 


Andrzej Miszczak 
Gdynia Maritime University 


ABSTRACT 


This paper presents an analytical solution of velocity components of unsymmetrical oil 
flow and pressure distribution in radial journal bearing gap for hydrodynamic unsteady 
lubrication with viscoelastic oil. Numerical calculations are performed in Mathcad 11 
Professional Program, with taking into account the method of finite differences. This me- 
thod satisfies stability conditions of numerical solutions of partial differential equations 

} and values of capacity forces occurring in cylindrical bearings. Exact calculations of pres- 
LA Wat! sure in journal bearing and its load capacity may be useful to prevent from premature 
wear tribological units of self ignition engines, especially those applied in ships. 


Key words : viscoelastic unsteady lubrication, analytical and numerical calculation, capacity forces, hydrodynamic pressure 


INTRODUCTION 


Correctness assessment of functioning the machines used for driving systems of various transport means e.g. diesel engines 
in which journal friction units are installed, depends to a large extent on assumed computational models, correct estimation of 
assumed appropriate simplifications, if necessary, and then on an assumed numerical method for the determining of operational 
parameters. During maneouvres of sea-going ships as well as when driving cars many frequent changes of engine operational 
parameters occur, especially of its rotational speed and loading. Similarly it happens during sailing the ship in rough weather 
when non-stationary loads on its propulsion engine happen. 

Both car vehicles and majority of sea-going ships are equipped with a propulsion system of a single self-ignition combustion 
engine. Seizure of trigological system of such engine is equivalent to depriving the vehicle or ship of its propulsion and thereby 
of its serviceability [3]. If such an event occurs during storm emergency situations or even a catastrophy may happen [4]. One of 
the ways to avoid such failures is to apply a lubricating oil of required properties, especially of appropriate viscosity and lubricity. 

Therefore the main aim of this work is to present analytical-numerical calculations of distributions of hydrodynamic pressure 
values occurring during non-stationary impulse lubrication of bearing surfaces at viscoelastic oil flow. Viscoelastic properties are 
characteristic for all oils which contain various bettering admixtures or in which some impurities such as lead salts, soot or dust, 
happen. All such impurities and admixtures are typical for land and sea transport. 

Taking into account the above mentioned observations it is necessary to precisely analyze the influence of non-stationary 
load impulses which are transferred through a propulsion system to a slide friction unit and result in such characteristic quantity 
of the bearing as its load-carrying capacity determined on the basis of pressure distributions. 


In Fig.1 the structure and loading scheme of a slide journal bearing is presented. 


i h(@,z,t) . 
GR 


s 
N 
LU 


i z i 
Fig. 1. Structure and loading scheme of a radial slide journal bearing and its characteristic dimensions 
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As shown in the scheme the bearing sleeve axis may undergo a skew respective to the journal axis. The quantity is described 
by the angle y. In such case the oil gap height depends on the tangential variable @ and the longitudinal coordinate z. In the case 
of non-stationary impulse loads the oil gap height additionally depends on time t. 


THEORETICAL MODEL 


Oil flow through the cylindrical gap height of radial slide bearing is described by the momentum conservation equations and 
continuity equation [1, 2, 6, 7, 8, 10, 11, 15]. Additionally, Rivlin-Ericksen constitutive relationships were assumed. The equ- 
ations in question are of the following form : 


Div S = pdv/dt , divv=0 , S= -pI +1 A; + a(A1) + BA2 (1) 
where : 
S - stress tensor p -oil density 
Div S - stress tensor divergence t -time 
v - velocity vector p -pressure 
divy - velocity vector divergence I - unit tensor 


A, and A) -two Rivlin-Ericksen strain tensors of three material constants No, o, B, 
where : 
No - dynamic viscosity ©, B - pseudo-viscosity constans of oil 
The coordinates of A4, Ag tensors are described by the symmetrical matrices defined as follows : 


A] = L +f L T 
A2 = grad a + (grada) +2L L (2) 
Lv + x 
a=LVT > 
ot 
where : 
L -tensor of oil velocity vector gradient s`! a - acceleration vector 
LT -tensor with matrix transpose s` grada -acceleration vector gradient 


The product of the Deborah and Strouhal numbers, marked DeStr, is assumedof the same order as the product 
of the Reynold's number, relative radial clearance and Strouhal number, marked RewStr. Moreover : DeStr >> De = aw/no . 
where : y - relative radial clearence, and w - angular speed of cylindrical bearing journal. 


The following is additionally assumed : 


rotational motion of the journal with the tangential speed U = wR 
unsymmetrical, non-stationary oil flow through bearing gap height 
non-stationary viscoelastic properties of oil 

constant oil density p 

the characteristic gap height h(@,z,t), in the cylindrical bearing 

no slip between bearing surfaces 

R - radius of cylindrical journal 

2b - length of the bearing in question. 


VVVVVVV V 


Neglecting the terms for the relative radial clearance y = €/R= 10° in the basic equations defined in the cylindrical coordi- 
nate frame : @, r, z, as well as taking into account the above mentioned assumptions, one can obtain : 


3 
Wo ___1 Sp, Me O/ Mo), BO Ve m 
ôt pRôp p orl ô ) pôt? 
0-2 (4) 

av 3 
Ov, _ Lop nm 2 ON, Bê (5) 
ôt pôz p orl ar) p atér? 

dv 

INe, N M2 2 (6) 


R op coro 
for:0<@<2m , -bSz<+b , OSrs<h , where: h - characteristic gap height. 


The symbols: Vo, Vr, Vz represent the respective oil velocity vector components : tangentially directed, that along gap height, and 
longitudinally directed. The following relationships between the dimensional and dimensionless quantities are assumed [12 , 13] : 


r= R(1+yrı) > z = bz; 2 t= toti ’ h= eh; ’ Vo = Yvo ’ Vr = UW vr 7 
vz = (U/L1) Vz1 5 P=PoP1 » Po = UNoR/e” , Lı =b/R o 
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Reynolds number, modified Reynolds number, Strouhal and Deborah number are assumed in the following forms : 
Re = pU£/no , Rey = POE Mo , Str=R/(Uto) , De = BU/MoR) (8) 
hence : DeStr = B/(Moto) = Des , RewStr = 087/(Note) = Res (8a) 


For the commonly applied inhibitors the coefficient B satisfies the inequality: 0 < B/ty < No. Values of the coefficient B vary 
from 0.000001 Pa's“ to 0.01 Pa's“. The dimensionless symbols are marked with the lower index "1". The equations (3) + (6) take 
the dimensionless form : 


3 
R = eS. Des — Ô Voi (9) 
ôt, Op Or, = P 
o=% (10) 
1 
3 
Rese oe circ Ni epee a "a (11) 
ôt, CZ, Ob \ Of, Ot, Or, 
Ov 
a na 1 Vy -0 (12) 
op oY Li OZ, 
for:0<S@<2n , -1<zı<+1 , 0<rı<hı 


GENERAL AND PARTICULAR SOLUTIONS 


A new variable is now introduced : 


D 
y=nN, Nat Bes » ty ors eng | (13) 
2%; A 


and solutions are assumed to have the form of the following convergent series [5] : 


Vol = Voor (XQ, Z; ) Van apa2] von (Xo Q, ae Foen (14) 
Va = Vzoz (X> 0521) + a Vas (X QZ) + [Pe Vz05 (Xo P 21) F aeee (15) 
Va = Vros 1,0.2)+ v Pz) [2]; Vias (X Qs Z1) ties (16) 
pı polozit) + py(qzt+[ 2) Pio (Ps Zis ty) ee (17) 


for:t;>0 , 0<Des<<1 , (Des/t})<1 


In the equations (9) + (11) the derivatives respective to the variables t1, rı are substituted for the derivatives relative to the 
variable x, by using the relationships given in App. 1. Therefore the variables tı, rı are substituted for the variable x. Now the 
series (14) + (17) are introduced to the equation system (9) + (11). Next, the terms multiplied by the parameters of the same 
power, (Des/t))§, are successively compared for k = 0, 1, 2, ... . Hence the following sequence of ordinary differential equations 
is obtained: 


2 

d Voor 5, Noor 1 Oro (18) 
dy? dy N? a 

Pyar oy Waz L Pio (19) 
dy? dy N?’ 62, 
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d?v d?v d?v 
T tE Avo = Pu r L a 
dy dy N^ 0 dy, 2 dy 
d?v is +2% dvs + 4V 5 1 Py dva ii dvs 
dy? ° ôz dy? 2” dy 
dvs 5) p22 gv _ 1 OP d’ Vou 4 t d Volz 
dy? dy p22 2 ap dy? 2 dy? 
dva +2% 725 FEV PEE Va 1 d “az 
dy, % N4 ôz dy, 2 dy 
and so on. 


The general solutions of the differential equations (18), (19) are of the following form : 


Vo0x(X0,Z1) = 


Vz2(%-Ħ.Zz1) = 


Co1V01(X) + Co2V02(X) + Voo3(X.,Z1) 
C21V01(X) + C22V02(X) + V203 (X-Ọ,z1) 


where: Cg) , Cg2 , Cz1 , C22 are integration constants. The particular solutions of the uniform 
and non-uniform ordinary differential equation are as follows : 


2% 2 
va) = eet » V(X) =l., eios dy, 
0 


V 903 (X.9,2Z,) = 2N? 


Vn Opi 


dr Di Q(x) 


Q >’ o> == 
(x) V203 (X-Ọ;Z1) ON? Oz, 


Op 


x 


x 
Ay) = few erfy dy, — erfy few dy, 
0 


For t;0 , one obtains N> œ , hence y> œ . If t; > œ , then N> 0, as well as for rı > 0 one obtains y> 0. 


If tı > 0 andr, = 


Voi(X) = 1°7/2 for: 


Voi(X) = 0 for 
Vio3(X) = 9 for 
Vo(X%) =0 for 


2 
-i Po 
2 ao 
2 
o 10 
We a 
2 OZ, 


The cylindrical journal executes only the rotational motion in the @ direction. Hence the oil velocity component 
on the journal surface in the tangential direction is equal to the velocity of the cylindrical surface of the journal. 


0 
where : 0 < %1 <S% = rN 
0 then x = 0. The following limits are true : 


y>% , t 90 , Noo 
7X0 , 1 =0 , 0<t<t <, N>0 


X0 , 1 =0 , O<t<th<o , N>O 5 where:i=@,z 


:X>0 , 1 >0 , t>” , N>0 


for: ¥-0 , 1, >0 , t>” , N-0 


for: Y-0 , r,>0 , t>» , N-0 


(20) 


(21) 


(22) 


(23) 


(24) 


(25) 


(26) 


(27) 


(28) 


(29) 


(30) 


The longitudinal component of oil velocity on the journal's cylindrical surface equals zero because the journal is motionless 
along the z-axis. However the shaft undergoes pulsatory changes of its location with time along the gap height direction. Hence 
the gap height is time dependent. Thus, the following boundary conditions appear : 

Voox(X =0)=1 , Veox(¥=0)=0 , vos(% = 0) = Str - (Ahj/dty) 
for: tr; =O GS Y=0 , 0<t<t <, N>0 


(31) 


The cylindrical sleeve surface is motionless in the tangential, longitudinal ( axial) and radial directions. Viscous oil flows 
around the sleeve. Hence the oil velocity on the sleeve surface equals zero in the tangential and longitudinal directions as well as 
in the gap height direction r. Thus the following boundary conditions are valid : 


Voox(X =M)=0 , veox(X=M)=0 , viox(x = M) = 0 
for: Th; & YONh, =M , 0<ti<t <% , N>0 
where : h = eh, - gap height , hı- dimensionless gap height 
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(32) 


Applying the conditions (31), (32) to the solutions (24), (25) one obtains : 


CoiVoi(X = 9) + Cor + Vonx(X=9)=1  for:rı=0 
CoiVoi(X = M) + Co2 Bg Vo03(X = M) =0 for: 1 = hı 
Cz1V01(X = 0) + C+ Vv203(% = 0) =0 for: 1; =O 
Cr1V01(X% = M) + Cy + V203(% =M) = 0 for: rı = hı 


(33) 


If the limits (29), (30) are taken into account the set of equations (33) yields the following solutions : 


1+ V9; (M M 
? 03 M) Cy = Yat) Cas1. Cy =0 (34) 
Vo, (M) Vol 


Now, to the right hand sides of the equations (20), (21) the solutions (24), (25), (26), (27), (28), (34) are substituted. 
Hence the general solution of the equations (20), (21) obtains the following form : 


Voix (X021) = Cos V11 (X) + Coa V12 (X) + Vora 52) (35) 
Vas (021) = C33 Y1 00) + Cra V0 %0) + Y 213 (X92) (36) 


where : Co3 , Coa, C23 , Cz4 - integration constants 
93 > @ 8 


The particular solutions are as follows : 


vu) = ye” (37) 


x 
Vin(X) = ye [ae qx) (38) 
5X1 


d2 
y2 1 
Vgi3 ZC) = sido fic ot X11 +2) — (l)e 7 iawn a Vir(K dx) + 


Pu 


tof [iste Sh ws tI 2 Ap 


; — Coin +2) pV dx 


(39) 
OP 


1 


d? 
Va3 ZC gt) = vof Caki +2)- (i)e areala Vio (Xi dx + 
1 


l OP 


N? &z, -= C aX 1 +2) Vi dx; 


X 
E vn | í + 1 Srl al 


for:0<ð<%ı SX 


The solutions (35), (36) represent corrections to the components of oil velocity because of its viscoelastic properties. 
On the basis of the solutions (37) + (39) for y> 0 , rı—> 0 , N > 0 the following limits are achieved : 


9 x p) 
an Vin (XY) = aim ye X [4 em Jen =-1 (40) 
l 


—>0,N>0 —>0,N>0 3 


The following limits are also true : 

vii(X) = 0 for: X0 , r =0,0<t<t <, N>0 
Vi(X%) =-l for: X0 , 1. =0,0<t<t <, N>0 (41) 

vi3(%) =0 for: ¥ 0 , 1 =0 , O<t,<th<0 , N>O ; wherei=Q,z 
The corrections to the oil velocity components cannot have the same boundary conditions as those previously assumed (31), 
(32) for the cylindrical journal and sleeve in the longitudinal and tangential directions. Therefore the corrections satisfy the 

following boundary conditions : 

Voiz(X% = 0) = 0 5 Vas(% = 0) =0 for: 7, =O SG y=0 „0<t <t <œ, N>0 


42 
Voix(X% = M) = 0 x Vzix(% =M)=0 for : rh; & Y-Nh; = M „0<t <t <œ, N>0 (42) 
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Applying the conditions (42) to the general solutions (35), (36) one obtains : 
CosvVi(X = 0) + Cosvia(X = 0) + Vois(X = 0) = 0 for: r; =0 
Cosvis(X% = M) + Coavi2(X% = M) + Vois(X =M)=0 for: 11 = hy 


43 
Ca3viil(% = 0) + CraV12(X% = 0) T V213(X% = 0) =0 for: 1, =0 we 
Cr3Vi(X = M) + Cravio(X = M) + vris(X=M)=0 for:rı =h; 
If the limits (41) are accounted for the set of equations (43) yields the following solutions : 
. =h, N 
Cg =- Virsa = hN) L ) , Ca =0 fo:i=Ọ,Zz (44) 


Vu (x = h,N) 
The general solutions of the oil velocity components (15), (14) at making use of the solutions (25), (35), (36), 
can be presented in the following form : 


Des Des 
Vol = Caio (x) + Coo + V 003 (%-9,Z,) + EN Cus 1 (x) T CoVi2 (x) + Vo13 (X-Ọ.Z; )| A of 28) (45) 
1 1 


Des Des 
Va = CaVoi 0 + Cz + V 203 (%0,z1) + Cav 1 + C32 O) + Ya Z )]+ (22) (46) 
1 1 


If tj ce then N->0 , hence % =rjN-0. 


For further analysis it is worthwhile to find the following limits : 


lim N? vin) = lim N° ye x <0 (47) 
y—>0,N>0 
x 
2 
lim N?v = lim N2ye* - |= er dy, =0 48 
x>0,N>0 1200 = ¥>0,N30 x : Ki (48) 
for:0<%ı<hıN , 0<ty<o , O<riı<hı , -l<z)<+l , 0<ọ<27 


VALUES OF OIL VELOCITY AND PRESSURE 
AT NON-STATIONARY NEWTONIAN LUBRICATION 


If viscoelastic properties of oil are neglected, then on the basis of the solutions (45), (46) 
the oil velocity components in the @ and z directions, at non-stationary flow, are of the following form : 


Jn Op erf(r, N) avn Op 
Voor (9.1 5Z1,t))=+1- pe ON oa Q(x = Nh,) FA 2N? w (x=Nr,) (49) 
Vn erf(r,N Jn 
Vros (Q.1 52151) = Pro ü= Nh, ) S. Pro olx T Nr,) (50) 


2N? az, erf(h,N) 2N? éz, 
fr:0<t <> , O<r<h, , -lSz) <1 , O<@<2m , O<S<X2SX%1SX%X=rnN<hN=M , hy =h,(Q,z)) 
The oil velocity components (49), (50) are now inserted to the continuity equation (12) 
and next the equation is integrated respective to the variable ry. 


The oil velocity component v,os in the gap height direction is not equal to zero on the cylindrical journal surface due to 
impulse displacements of the shaft. Therefore by applying the condition v,oy = Strdh,/ot, for rı = 0, the following form of this oil 
velocity component is obtained : 


"Ov "0 Vaos dh, 
309 ° Le du At, 
for:0<t;<o , O<m<ry<sh, , -l<z)<+l , O<@<2m , O0<S<X2<%1<X=rıN <hıN=M 

The oil velocity component v,ox equals zero on the sleeve surface. Integrating the continuity equation (12) along the direc- 
tion r and applying the boundary condition (32) for rı = hy, to the oil velocity component in the gap height direction, and making 
use of the conditions (29) one obtains the equation : 

a" h, 


1 ə oh 
— |v alt, +—— Iv adr, = Str — 
z posh E T J zoz Gry at, (52) 
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The solutions (49) + (50) are now inserted to the equation (52). One then obtains the following modified Reynolds equation : 
h; 
erf (rN )dr ; 
fa || SON ‘ Pio 
Bue Olx =Nh,)— [Q(z = Nr, ar, + 
2N? G0 || erf(h,N) Oy 
h; 
erf (r N)dr 
m ola ; AP 
+ — Q(x = Nh,)- [om = Nr, )dr, = (53) 
2N°L? ôz, || erf(h,N) > az, 


hı 
=- 1 SfN) dr + Str oe 
erf(h,N) Ot! 


for:0<m<r<hy , OS@<2n , -1<z1<+1 , 0<ti <% , OSXISXSUN , 0S N(t}) =0.5(Res/t})?°< 00 


The modified Reynolds equation (53) defines an unknown pressure function pj (Ọ, Z1, t1). If the dimensionless time tı 
approaches infinity, i.e. the coefficient N approaches zero, then the equation (53) approaches the classical Reynolds equation 


(see App. 2): 
EA h? OP 10 +12 h Pro 6 F954 (54) 
a dp ) Li ôz Oz, a) at, 
for:0<@<2n , -l<z,<+l 


The dimensionless time-dependent height of the bearing gap height, at accounting for its periodical disturbances, 
may be described by the following relationship : 


hy =[1+A-cosp+s-z-cos(@)] - f(t1) ; fti) =[1 + c: exp(-ty ` ty - Mo)] (55) 
where : s= “tan: skew coefficient. 


The symbol œo stands for an angular velocity given in [s1] which describes the disturbances in the gap height direction for 
unsteady oil flow through the cylindrical bearing gap height, and "c" stands for a coefficient used to control values of gap height 
changes. If c- value is positive the bearing gap height is increased, if negative - the bearing gap height is decreased relative to the 
classical gap height. If tı approaches infinity then the gap height equation (55) approaches the classical gap height equation 
independent on time during stationary motion. 


PRESSURE CORRECTIONS FOR VISCOELASTIC OIL PROPERTIES 


Particular solutions of oil velocity components in the directions @ and zı due to viscoelastic properties in non-stationary 
motion are contained in the solutions (14) , (16) multiplied by the term DeStr/t;. By making use of (38), (39), (44) and the 
boundary conditions (42), the corrections of the oil velocity components (35) , (36) can be transformed to the following forms : 


hN 
Des 4 Ry? | O f 
PS y a (@sZisisty) = —oeNre pnt) Pu [x2 (x)dx + 
ti pe ð qN 
h,N hN 


2 2 
OOd -Q x =hN) f We ™ Q, (dy + 
N 0 


N f2 2 OP 10 
+~ (k -12)0,(¢=h,N) |+ 22 
5 ( 1 2) ıx=h; | a 


q 


rN h,N 
1 2 p ÔPi0 1 
+Q =N) [xie ™ Q, dy | -———— OOG + (56) 
=à [r 200K TN Go i (23 Ody 
hN nN 
-Q(x =h N) JQ3(WQG)dx + 24H = 1N) [202x )dx |f + 
0 0 = 
BBN nein An o i 
-= _ Ia (y=, N) [Q3(y)dy- 2, (y= 4,N) odr- aOd 
Jive’ eeu) NK =A; fosoa ı X= Joswir- foa 3\X)0X 
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Des 4B 1 T 
Dea et a 
ti pe Li OZ, rN 


N? ôp h,N 
-roaa å z “la, (q)dx - Q(x = hy N) Jue CO, (Ody + 
Z] qN 
No, ; P biN (57) 
+Q =N) Jue” woal- E [2123 @QQ)dx+ 
1 qN 
h,N 
-Q(x =h,N) fo (1) Q(x)dx +Q (x = WN) Jo; a 
where 
t] 2 2 aà Xo a 
2,(y)= lee dy, , TEE Ejen " , wha x +2) (59) 
8 Al 0 


and:0<t;<o , OSry<r, Sh, , -lSz,<1 , O<@<2n , 0<% <SX%X=rıN<hıN=M 


The corrections (56) , (57) are now inserted to the continuity equation (12) 
and then both its sides are integrated respective to the variable r. 


From the viscoelastic oil properties it results that the corrections of the oil velocity component along the gap height must equal 
zero on the journal surface at rı = 0. Hence the correction of the oil velocity component along the gap height is as follows : 


e tai 
Varl 2.4, h) =— [Vor (9, 2,,7,t,)dy |+——— [van(@. Z 51), t,)dr, (59) 
OP | 6 Li 62 | à 

The velocity vector corrections cannot change the boundary conditions (31) , (32) which are assumed on the journal and 
sleeve surfaces in the direction of the bearing gap height. Therefore the oil velocity vector corrections in this direction are equal 
to zero on the sleeve surface at r;=h,. By applying this condition to the solution (59) the modified Reynolds equation can be 
obtained : 


hı m h,N 272 
22u pe f Kak aoka) a +(e ¥8 )n?0,6=h) = 


tN 


1 a la hy hN 
t= ne^ fho ) da-ra =N ) dr + th etn’ )n20,(,=hN) = 


aT N 


2N h, : h,N h,N 
=— 4 fret o (x =N) fa. (wdx - Q(x = 4N) Jo, (x)dx - Jo, ()Q, er lins 


erf (Nh, ) ; nN 
hN 
p O Pio ign i (60) 
p re 2 (x) Qa(y) Q(y) dy + 
Jmerf(h,N) 20 | dp j J a ee 


nN = -h?N? hy 
+O(L=HN) [AWOBO a dr = N =hN) [OOB a + 
0 0 


2 1 8 ap h; hı N 
-o 10} fren QOO 230) 20) dy + 
"MR erf(h,N) L? 6z, | 6z; J ' ee Rene 


rN -h? N? h; N 
+Q (x =N) Jo 0) Qs cnn ea (x = h,N) jo 09204) dx |+ 
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a | ap hy hN h,N i 
aye] fen] JQaGcdx +2,(2=HN) fo)e dx |an + 
p p 0 rN 0 
jeter u 2 
-— z7 ale= BN) [Orne * dx Jr + (60) 
0 
h h,N h,N 
1 8 | apo) e" p 
ERA az Ire į Jew) dy + Q,(y =N) fot) * dy | dr, + 
ie n 2 
ye Bl = bw) 2.x) xe* dy 
0 
for:0<m<rsh; , 0<Q<2n , -1<z<+1l , O0St)<0 , 0S ¥)<XSh|N , 0<N(t) =0.5(Res/t1)°> < o 


The modified Reynolds equation (60) determines unknown functions p11(Ọ, Z1, t1) of the corrections of pressure values, 
resulting from viscoelastic oil properties during non-stationary motion. 


NUMERICAL CALCULATIONS 


The dimensionless distributions of values of the pressure pọ and its corrections p11, p12, ... in the lubrication area are determi- 
ned by means of the Reynolds equations (53) , (60) with using the gap height (55). On the boundary of the area, dimensional values 
of the pressure and its corrections obtain the value of the atmospheric pressure pat. The lubrication area and the gap height are 
shown in Fig.2. In order to determine a dimensional value of the gap height the dimensionless gap height indicated in Fig.2 should 
be multiply by the radial clearance £. The lubrication area is defined by the following inequalities: 0<@<am , -l<z, <1. 


hy hy 


9 


Fig. 2. The example lubrication area with the time-variable height of lubrication bearing gap 


The values of the dimensionless time tı increasing from 1 to 2000 and even more express the time of departing from the 
instant of impulse occurrence of the force acting on friction unit. The impulse which causes the gap height decreasing was 
denoted with the negative value of the gap height change coefficient c = -1/4, and that causing the gap height increasing — with 
the positive value of that coefficient c = 1/4. 


In order to calculate definite hydrodynamic pressure values the following input data were assumed : 


¢ relative radial clearance : y =10° ¢ shaft radius : R =0.08 [m] 

¢ oil viscosity : No = 0.03 [Pas] + dimensionless bearing length : L,=1 

è il pseudo-viscosity coefficient: B = 6-107 [Pa:s?] + skewness coefficient : s =0.05 

¢ oil density : p = 900 [kg/m] + angular speed of sleeve perturbation: w,= 0.2 [s!] 

¢ shaft angular speed : n = 1500 [rpm] ¢ characteristic time : tə = 0.001 [s ] 

¢ relative excentricity : à =0.5 + time intervals : tı =1 ; 100; 10000 ; œœ 


For the so assumed data and on the basis of the Reynolds equation (53) the dimensionless values of pressure distribution were 
numerically determined for the gap height defined by the equation (55) with the use of the method of finite differences and the 
software Mathcad 11 ( Fig. 3, 4, 5, 6 ). In order to obtain real dimensional distributions of pressure values, the calculated dimen- 
sionless values of pressure distributions should be multiplied by the dimensional coefficient URy,/€°. The distributions of dimen- 
sionless hydrodynamic pressure values were presented for the dimensionless time intervals t = 1; t; = 100; t; = 10000; tı = œ and 
the bearing gap height change coefficient c =+ 1/4. 
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Fig. 3. Distributions of the dimensionless hydrodynamic pressure in the cylindrical bearing gap for the dimensionless time t]=1 counted from the impulse 
instant and for the gap height change coefficient c = +1/4 (lubrication gap height increased due to impulse load), and for c = -1/4 ( lubrication gap height 
decreased due to impulse load) at accounting for a skew of the journal 


(fe 
ELE 


fi 
Aa. 
re 
ji ieee 
Fi 


rar, 
aa, 


"a 
ete: 


U7 


0.0 


Mwy, 
(7 
pre 
7 


SARS 


PT = pio Omi W 

3.0 E NN ‘ 
Lj=1 |! 
à =0.5 i 

204 s =0.05 f | 
ca P10max = 1.227 fi \ { i 


j 
Í 
m 


SSS 
SS 
aa 
LS 


ts 
‘ou, 


2 
CE roe ee 
SS 


2n/3 a Ọ 


Fig. 4. Distributions of the dimensionless hydrodynamic pressure in the cylindrical bearing gap for the dimensionless time t;=100 counted from the impulse 
instant and for the gap height change coefficient c = +1/4 (lubrication gap height increased due to impulse load), and for c = -1/4 (lubrication gap height 
decreased due to impulse load) at accounting for a skew of the journal 
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Fig. 5. Distributions of the dimensionless hydrodynamic pressure in the cylindrical bearing gap for the dimensionless time t;=10000 counted from the im- 
pulse instant and for the gap height change coefficient c = +1/4 (lubrication gap height increased due to impulse load), and for c = -1/4 (lubrication gap 
height decreased due to impulse load) at accounting for a skew of the journal 
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Fig. 6. Distributions of the dimensionless hydrodynamic pressure in the cylindrical bearing gap for the dimensionless time t; = œ counted from the impulse 
instant and for the gap height change coefficient c = +1/4 (lubrication gap height increased due to impulse load), and for c = -1/4 (lubrication gap height 


decreased due to impulse load) at accounting for a skew of the journal 
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For the same data on the basis of the equation (60) the numerical calculations of the dimensionless values of hydrodynamic 
pressure corrections which result from oil viscoelastic properties, were performed. Their results are presented in Fig. 7, 8, 9 for 
the dimensionless time intervals tı = 1, 2, 10. For tı = 100, tı = 10000 and tı = œ, the calculations of dimensionless values of the 
hydrodynamic pressure corrections were also performed, but they have not been attached here as being negligible. To obtain 
dimensional values of the pressure corrections the dimensionless values shown in Fig. 7, 8, 9, should be multiplied by the 


dimensional coefficient URy,/€°. 
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Fig.7. Distributions of the dimensionless values of the hydrodynamic pressure corrections resulting from viscoelastic properties of oil in the cylindrical bear- 
ing gap for the dimensionless time t; = 1 counted from the impulse occurrence instant and for the gap height change coefficient c = +1/4 (lubrication gap 
height increased due to impulse load), and for c = -1/4 (lubrication gap height decreased due to impulse load) at accounting for a skew of the journal 
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Fig.8. Distributions of the dimensionless values of the hydrodynamic pressure corrections resulting from viscoelastic properties of oil in the cylindrical bear- 
ing gap for the dimensionless time t; = 2 counted from the impulse occurrence instant and for the gap height change coefficient c = +1/4 (lubrication gap 
height increased due to impulse load), and for c = -1/4 (lubrication gap height decreased due to impulse load) at accounting for a skew of the journal 
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Fig.9. Distributions of the dimensionless values of the hydrodynamic pressure corrections resulting from viscoelastic properties of oil in the cylindrical bear- 
ing gap for the dimensionless time t; = 10 counted from the impulse occurrence instant and for the gap height change coefficient c = +1/4 (lubrication gap 
height increased due to impulse load), and for c = -1/4 (lubrication gap height decreased due to impulse load) at accounting for a skew of the journal 


In Fig. 6 one can observe that when a load impulse occurs sufficiently far in time from the impulse occurrence instant, i.e. 
when tı—> œ, then the distributions of pressure values approach the pressure distribution identical as regards its values and shape, 
both at the impulse increasing the gap height c > 0 and that decreasing the gap height c < 0, which can be also achieved from the 


classical Reynolds equation (54). 


From the analysis of the pressure corrections (due to viscoelastic oil properties) it results that only in the initial instant tı = 1 
after impulse occurrence (Fig.7) the corrections really influence the total pressure value. At the so assumed time instant t; = 1 the 
share of the corrections of pressure p11 in the value of the basic pressure pjg amounts to about 6% (at the gap height decreased 
due to impulse load, c = -1/4) and to about 18% (at that gap height increased due to impulse load, c = +1/4). The values were 


calculated for the relevant maximum values shown in Fig. 3 and 7. 
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DISCUSSION OF RESULTS AND CONCLUSIONS 


O Numerical analysis of hydrodynamic pressure values during the unsteady impulse loading of radial slide cylindrical bearings 
reveals that in the first time interval counted from the instant of impulse load occurrence very high changes of the hydrody- 
namic pressure may appear, and also very high changes of bearing loads in comparison with their load carrying capacities 
which shall occur at no impulse load. 

O In the case if due to an impulse load the journal changes its location relative to the sleeve in such a way that the lubrication gap 
height increases (left column of Fig. 3, 4, 5 and 6) then the bearing will suffer sudden drop of its load carrying capacity by 
a few dozen percent. As time runs after the impulse load occurrence the hydrodynamic pressure in the bearing gap height 
increases up to its value occurring under regular load (without any impulse). 

O In the case when an impulse load results in decreasing the gap height (right column of Fig.3, 4, 5, 6) then an increase of 
pressure values appears in the initial phase of impulse loading and next, as time runs, the pressure decreases down to the 
hydrodynamic pressure value relevant for the bearing under regular load (without any impulse). 

O The mixed case may also happen when the journal displaces itself due to simultaneous occurrence of two impulses of 
opposite tendencies leading to the decreasing and increasing of the gap height relative to its initial location. Then, rises and 
drops of hydrodynamic pressure in comparison to its initial value, may happen. Such hydrodynamic pressure changes may 
lead to an accelerated wear of elements of the cylindrical slide friction units in question. 

O The accounting for the impulse-load-induced pressure changes in designing the cylindrical slide friction units, would contri- 
bute to elimination of engine failures resulting from seizure of cylindrical slide bearings in the service conditions in which 
impulse loads often occur. Transport safety would be this way improved. 
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The equation (53) at N-0 approaches the following form : 
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After realisation of the calculations the classical Reynolds equation is obtained in the form of (54) for a cylindrical system. 
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NOMENCLATURE 

a - acceleration vector [m/s2] 

Ay - Rivlin-Ericksen strain tensors [s-!] 

A> - Rivlin-Ericksen strain tensors [s-2] 

b - a half of bearing length [m] 

c - coefficient for controllling gap height changes 
Cr, C22, C33 , Cra, Cots Cor, Cos, Cong - integration constants 
De - Deborah number 

h - gap height in the cylindrical bearing [m] 

hı - dimensionless gap height 

I - unit tensor 

L - tensor of oil velocity vector gradient [s-!] 

LT - tensor with matrix transpose [s-!] 

L, - dimensionless bearing length 

N - dimensionless number 


- estimate of all remaining corrections of velocity 
and pressure components 


p - pressure [Pa] 

Po - characteristic value of hydrodynamic pressure [Pa] 

Pi - total dimensionless hydrodynamic pressure 

P10» P11; P2  - dimensionless corrections of hydrodynamic 
pressure 

P - load [N] 

T - radial coordinate [m] 

Tı, f2 - dimensionless radial coordinate 

R - radius of cylindrical journal [m] 

Re - Reynold's number 

s - skew coefficient 

S - stress tensor [Pa] 

Str - Strouhal number 

t - time [s] 

fo - characteristic time [s] 

ti, to - dimensionless time 

U - tangential journal velocity [m/s] 


v - velocity vector [m/s] 


Vgs Vra Vz - dimensional values of tangential, radial and axial 
components of velocity vector, respectively [m/s] 
Vọ1 s Vri > Vzi - dimensionless values of tangential, radial and 


axial components of velocity vector, respectively 

Vox» Vroz» Vzox - dimensionless components of oil velocity vector, 
without accounting for changes due to disturbing 
impulse 

Vold> Vrid.> Vzid» V22; V22» Vzax- dimensionless corrections of oil 
velocity vector components, resulting from 
disturbing impulse impact on a bearing at 
sufficiently close instant from the impulse 
occurrence 


Vio» V20 - parts of dimensionless velocity vector 
components, dependent on shaft rotation, without 
accounting for disturbing impulse action 

V11, V12 - parts of dimensionless velocity vector 
components, dependent on shaft rotation, with 
accounting for disturbing impulse action 

Vi903 » Vz03 - parts of dimensionless velocity vector 
components, resulting from pressure gradient 
influence, without accounting for disturbing 
impulse action 

Vo13» Vz13 - parts of dimensionless velocity vector components, 
resulting from pressure gradient influence, with 
accounting for disturbing impulse action 

Z - longitudinal coordinate [m] 

Z] - dimensionless longitudinal coordinate 

(o APA e) - pseudo-viscosity constants of oil [Pa-s2] 

Y - skew angle 


- a value close zero 

- radial clearance 

- characteristic value of oil dynamic viscosity [Pas] 
- relative eccentricity 

- oil density [kg/m3] 

- tangential coordinate 

- dimensionless coordinates 

- relative radial clearance 

- angular journal velocity [s-!] 

- angular speed of sleeve perturbation [s-!] 
Q, Qi, Q,, Q, - auxiliary functions 
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Seating of machines and devices on foundation 
chocks cast of EPY resin compound 


A valuable book written by Prof. Karol Grudzinski and 
Dr Eng Wiesław Jaroszewicz, has been recently published. 


The book presents a modern method for the seating of 
marine and land-based machines and devices on chocks cast 
of EPY resin compound specially developed for this purpo- 
se. General requirements referring to the seating of machi- 
nery on foundations (especially those used in shipbuilding) 
are listed together with relevant evaluation criteria. The pro- 
perties of resin compounds used for foundation chocks, the 
background of chocking arrangement design and the techni- 
ques used for casting the chocks in place are also outlined. 


Many examples of so installed machines and devices are 
described, illustrating various possible applications of EPY 
compound to the seating of new machinery and the repairs 
of existing one. The results and descriptions of research aimed 
at finding solutions for many practical problems in this field, 
constituting a scientific basis of the methods developed for 
the seating of machines on their foundations, are also given. 


The book of 186 pages consists 
of the six topical chapters : 


1. Characteristics and types of chocking arrangements used 
for ship machinery 

2. The resin compounds used for ship machinery founda- 
tion chocks 

3. Design of machinery chocking arrangements with EPY 
compound chocks 

4. The technology of machinery seating on EPY compo- 
und chocks 

5. Applications of EPY compound chocks for the seating 
of machinery - practical examples 

6. Research on resin compounds used for foundation chocks 
(this chapter occupies amost a half of the book and is 
illustrated with 74 fotoghaphies and drawings) 


The attached bibliography contains 83 references. The 
book is ended with "The chronological list of research re- 
ports concerning Polish resin compounds used for founda- 
tion chocks, and their practical application for the seating of 
machinery"(containing 132 items both published and unpu- 
blished). 


The book is addressed to designers and shipbuilding 
technology specialists as well as the engineers and techni- 
cians dealing with the design, modernisation and execution 
of various heavy machinery installations on land. It may 
also be of use for the scientific workers and students of the 
university faculties engaged in the fields of shipbuilding 
and offshore technology, machinery design and maintenan- 
ce, industrial constructions and the building of roads and 
bridges. 


The first of the authors of the book, Prof. Karol Grudzin- 
ski, has worked at the Chair of Mechanics and Principles of 
Mechanical Design, Faculty of Mechanical Engineering, 
Technical University of Szczecin. The area of his scientific 
interests covers the mechanics of machines, mechanics of 
contact between real surfaces and tribology, with main fo- 
cus on experimental research, modelling and calculations of 
constructional joints. The results of his research has been 
presented in this book. Prof. Karol Grudzinski is an author 
and co-author of about 200 scientific papers. 


Dr Eng Wiesław Jaroszewicz, the other author of the 
book, has been angaged in the machinery seating technolo- 
gy since the early 1980s. He has owned the Marine Service 
Jaroszewicz company since 1990, which has developed the 
EPY resin compound and carried out machinery seating ope- 
rations based on its use. He is a co-author of many scientific 
publications, patents and reports. Over 7500 shipboard and 
land machines have been installed in the years of 1974 - 
- 2002 under his supervision and with his cooperation, in- 
cluding over 1500 main engines of ships. The practical exper- 
tise gained during this multi-year activity has been compre- 
hensively presented in the book in question. 


This English edition of the book 
is a translation from the Polish edition of 2002. 
The book assigned with ISBN 83-87377-94-5 mark, has 
been published by ZAPOL Publishing Company, Szczecin. 


KAROL GRUDZINSKI, WIESLAW JAROSZEWICZ 


POSADAWIANIE MASZYN | URZĄDZEŃ 
NA PODKLADKACH FUNDAMENTOWYCH 
ODLEWANYCH Z TWORZYWA EPY 


SZCZECIN 2002 


22 


POLISH MARITIME RESEARCH, No 2/2004 


MARINE ENGINEERING 


Fatigue "safe-life" criterion for metal elements 
under multiaxial constant and periodic loads 


Janusz Kolenda 
Gdańsk University of Technology 
Naval University of Gdynia 


ABSTRACT 


Periodic stress with Cartesian components given in the form of Fourier series is conside- 
red. To account for the mean stress effect the generalised Soderberg criterion for ductile 
materials is employed. An equivalent stress with synchronous components is defined by 
means of the equivalence conditions based on the average strain energy of distortion. The 
fatigue "safe-life" design criterion is formulated which covers the conditions of both static 
strength and fatigue safety and includes material constants that have simple physical in- 
terpretation, can be determined by uniaxial tests, are related directly to the applied loads 
and can reflect material anisotropy. 


Key words : design criteria, multiaxial loading, periodic stress, mean stress effect 


INTRODUCTION 


Machinery parts and structural elements are frequently sub- 
jected to simultaneous action of constant and cyclic stresses. 
In marine floating objects, constant stresses are mainly caused 
by deadweight and hydrostatic pressure, whereas those enco- 
untered in machinery systems are produced by torque, thrust, 
centrifugal force, etc. In the present paper also cyclic stresses 
induced by time-varying loads of periodic character are consi- 
dered. 

Each Cartesian component of the cyclic stress can be cha- 
racterised by its mean value and parameters of its alternating 
part. It is important to remember that the total strength of an 
engineering element is altered if residual stresses are present. 
Since residual stresses have a similar influence on the fatigue 
behaviour of materials as that of mechanically imposed con- 
stant stresses of the same magnitude [1], no distinction will be 
made between constant, mean and residual stresses. 

In order to ensure a fatigue safe life (theoretically infinite) 
under bending stress of mean value cp and amplitude ap two 
approaches can be indicated [2]. The first is based on the con- 
dition : 


a, =F, (1) 
where : 
ay =a TWh Cph (2) 


is the amplitude of the equivalent zero mean stress, 


Wp =0.1+0.2 -the asymmetry sensitivity index at bending 
Fy - the fatigue limit under fully reversed bending. 


Analogous conditions can be applied in the case 
of axial force or torsional moment. 


In the second approach the effect of mean stress is descri- 
bed by a "failure diagram" or by one of the empirical equations, 
such as Goodman's, Gerber's, Soderberg's or Bagci's equation, 
depending on a given situation [2+5]. In the following the So- 
derberg's equation for ductile materials is used in the form : 


oe (3) 


where : 

B - maximum allowable stress amplitude in fatigue 
"safe-life" design under asymmetric bending 

Re - tensile yield strength. 


BACKGROUND 


Introducing the safety factor : 


-B (4) 
ay 


and partial safety factors : 
R F, 


e a= 
ay 


S Cy 


(5) 
one obtains from (3) : 


f=£,(1-£,') 6) 
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So, the design criterion reads : 


f>1 (7) 
L.e., 
f-'+f;'<1 (8) 


Here the subscript "s" stands for static 
and "d" for dynamic parts of the applied stress. 


Obviously, the condition (1) concerns exclusively fatigue 
endurance of the material subjected to combined constant and 
cyclic loads, whereas satisfaction of (8) guarantees that not only 
the static strength of the material is not exceeded but also that 
the combination of constant and cyclic loads will not lead to 
fatigue failure. Another essential point is that the partial safety 
factors in (8) can be analysed and/or influenced separately. 
Therefore (1) and (2) is not used hereunder. 


The explicit forms of (8) read : 


Cc a 
—b 4b <] (9) 
R e F b 
for asymmetric bending, 


Ka + da <1 
e a 
for asymmetric push-pull force, and 


c a 
—t +-t<] 
R OR 
for asymmetric torsional moment, 
where : 


R, =R, 


3 


F, -1r 


3 


applicable for steels [4]. 


(10) 


(11) 


(12) 


(13) 


Here the subscripts "a", "b" and "t" denote axial, bending 
and torsional load cases, respectively, and F, is the fatigue li- 
mit under symmetric tension-compression. 

The aim of this paper is to extend the use of eqation (8) to 
multiaxial non-zero mean periodic stresses. For this purpose it 
is necessary to determine a reduced stress, equivalent in terms 
of fatigue performance of the material under the multiaxial 
stress. For example, in the case of in-phase stress with Carte- 
sian components : 


o,(t)=c,+a,sinwt i=x,y,z,xy,yz,zx (14) 


such reduced stress can be calculated for ductile metals 
by means of the average-distortion-energy 
strength hypothesis [6] as : 


Gq (t) = Coq + eq sinat (15) 
where its mean value and amplitude are given by : 
1/2 
cy +c, +e, — CC — CC, + 
Coq 2 2 2 (10) 
—C,C, +3(c2, eC c2) 
1/2 
az +a, +a, -—a,ā,-a,a, + 
aeq” (17) 
q 
—a,a, + 3(a2, + ai, + a2 ) 
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For the sake of brevity the stress components 6;(t), Oyz(t) 
and 67;(t) are dropped. It is noteworthy that in the case of out- 
-of-phase stress components : 


o,(t)=c¢, +a; sin(ot +B, ) 1=X,y,xy (14a) 
the average-distortion-energy strength hypothesis yields : 
1/2 
Aeg= la? +a, aya cos(B,—B,) +3a2, | (17a) 


According to the aforementioned hypothesis, the reduced 
stress (15) corresponds to that in an element of the specimen 
under uniaxial tension-compression test. Consequently, (10) 
can be used, which leads to the partial safety factors : 


R F 
fee f= (18) 
ST Cay d aeq 
and to the fatigue "safe-life" criterion : 
1 1/2 
—(c? + e = CC +3c2, ) + 
R, (19) 
I 1/2 
+h? tay =a,a, +3a2, <1 
i 
With (12), (16) and (18) one gets : 
-1/2 
l 2 2 2 
f, = as + Cy -c c, J+ REOS (20) 


Since normal stress components can be produced by loads 
of different mode (tension, compression, bending), and shear 
stress components — by torsion or shear, and the material may 
exhibit certain degree of anisotropy, the following modifica- 
tion of (20) is suggested [7] : 


-1/2 
R ÝR y 
Sc, |+| — Gs | + 
1I R, ge 
=a + 
Re R. R, 
f, = R “xR Cy (21) 
x y 
2 
Rete ~ 


Here Rx is the yield strength under static load relevant to 
the mean stress component cx. The remaining material con- 
stants are defined analogously. Equation (21) gives : 


-1/2 

c CC l 

f= 5 , 4 1=X,y,Xy (22) 
i R; R,R, 
With (13) and (17), the equation (23), similar to (22), 
can be written for the partial safety factor fy [8] : 
2 -1/2 
a,a 

a: x R 

f= aje i=xX,y,xy (23) 
: py FE] RF, 


where F; is the fatigue limit under fully reversed load relevant 
to the stress amplitude aj. Equations (8), (22) and (23) yield the 
following criterion of fatigue "safe life" under combined mul- 
tiaxial constant and in-phase loads [9] : 


1/2 
S C C,Cy P 
i R; R,Ry 
1/2 (24) 
+| > je ej 
i F, FE, 


Equation (24) may be called the generalised Soderberg cri- 
terion of an infinite fatigue life under non-zero mean in-phase 
stress. 


EQUIVALENT STRESS UNDER 
MULTIAXIAL CONSTANT 
AND PERIODIC LOADS 


Let us consider the stress with Cartesian components 
given by Fourier series : 


` : . (25) 
o;(t) =C; + Yay sin (pat + Bip) 1=xX,y,Xy 
p=l 
where : 
Cj - mean value of i-th stress component 


- amplitude and phase angle of p-th term in 
Fourier expansion of i-th stress component 

(Wo = 27/To - fundamental circular frequency 

To - common period of the stress components. 


dip» Bip 


Guided by the above presented criterion for the stress com- 
ponents (14), we shall try to model the stress components (25) 
by the equivalent stress components : 


(26) 
of) (t) = ole) + (ea) sin lot +0; ) i=x,y,xy 
where : 
cfa) - mean value of i-th equivalent stress component 
ala), @, - amplitude and phase angle of i-th equivalent 
stress component 
Weq - equivalent circular frequency. 


For determination of parameters of the equivalent stress 
components the theory of energy transformation systems [10] 
is used. According to this theory, a reduced stress model and 
a given multiaxial stress can be regarded as equivalent in terms 
of fatigue life of the material if during the service life the inter- 
nally and externally dissipated energies per unit volume in the- 
se two stress states are respectively equal. Under the assump- 
tion that the externally dissipated energy is proportional to the 
average strain energy of distortion, the following equivalence 
condition can be written [11] : 


= J. (t)at a folat 27) 
where : p ( i 
(eq) (eq) 
l+v| 0 + \O + 
dat) ” 


— olga) 4 309)" 


X 


is the strain energy of distortion per unit volume 
in the equivalent stress state, and : 


10- +05 -00y +302 5 


is that in the actual stress state, and : 
E - Young modulus 
v -Poisson's ratio. 


If the condition (27) is fulfilled 
then also the following equations are satisfied : 


Ty 2 Ts 
— SBO] a= fo Oa 
0 T, 0 


(28) 
T, 
ie 
a | o&)(t)o(t)dt = 
i (29) 
=— [o, (to, (t)dt 
o 0 
When : Mg = kop (30) 
where : k - a natural number, is assumed 
then (25) through (29) yield : 
ole =< (31) 
P 1/2 
eq) _ 2 
és ( q) 5 a, mA 
p=! 
alda (ea) cos (o, ~0, )= 
(33) 


-$a pay £08 (Byy~Byp) 


Of course, in design for an infinite fatigue life 
the evaluation of equivalent circular frequency can be avoided. 


FATIGUE "SAFE-LIFE" CRITERION 


Having determined the parameters of equivalent stress com- 
ponents, one can make use of the criterion (24) which, in parti- 
cular, is suitable also for non-zero mean stress with synchrono- 
us components of the form (26). In this instance one obtains : 


(ea) : cedola 
LIR) 


R,R, 
3 CAC 1/2 
$ y alo adad cos(p,- o, ) <] 
EE, 7 


(34) 


i i 


so that the fatigue "safe-life" design criterion for metal ele- 

ments under multiaxial constant and periodic loads becomes : 
1/2 

ec, 


RR, (35) 


1/2 


Yayay cos(B,, ~ By») <1 


xFy pal 


dzie 


i p=l 


Its extension to three-dimensional cases is straightforward. 
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CONCLUSIONS 


O The fatigue "safe-life" design criterion covering the condi- 
tions of both static strength and fatigue safety of metal ele- 
ments under multiaxial constant and periodic loads, was 
formulated. 

O The presented criterion includes material constants which : 
have simple physical interpretation, can be determined by 
uniaxial tests, are directly related to the applied loads, and 
can reflect the material anisotropy. 


NOMENCLATURE 
aa, ap, a - Stress amplitude under asymmetric axial force, 
: bending moment and torsional load, respectively 
ap - amplitude of the equivalent zero mean stress under 
asymmetric bending 


aj, Ci - amplitude and mean value of i-th stress component 
(i =X, y, Z, XY, YZ, ZX) 
aeq» Ceg  - amplitude and mean value of the equivalent stress, 


respectively 
aed), ofa) - amplitude and mean value of i-th equivalent stress 
component, respectively 


ain - amplitude of p-th harmonic in Fourier expansion 
of i-th stress component 
B - maximum allowable stress amplitude in design 


for an infinite fatigue life under asymmetric bending 
Ca, Cp, CŒ - mean stress value under axial force, bending moment 
and torsional load, respectively 


E - Young modulus 
f - safety factor 
fy, f - partial safety factors 


F,, Fp, F, - fatigue limit under fully reversed tension-compression, 
bending and torsion, respectively 
Fi - fatigue limit under fully reversed load relevant 
to the stress amplitude a; 
- natural number 
- tensile yield strength 
- yield strength relevant to the mean stress component c; 
- shear yield strength 
- time 
- stress period 
- phase angle of i-th component of the out-of-phase 
stress 
- phase angle of p-th harmonic in Fourier expansion 
of i-th stress component 
- Poisson's ratio 
- i-th stress component 
- equivalent stress 


Boda 


Gs 


© 


P PAT 
gri 


Aa aAas< 
Ag 
D 


- i-th equivalent stress component 

- phase angle of i-th equivalent stress component 

- strain energy of distortion per unit volume in the actual 
and equivalent stress states, respectively 

- asymmetry sensitivity index in bending 

- circular frequency 

(ON - fundamental circular frequency of the actual stress 

- equivalent circular frequency 


$ 
Poa 
o 
Ko} 


SS 
o 
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a C onference - 


Maritime Traffic Engineering 


On 20 and 21 November 2003, 10th International Con- 
ference devoted to problems of maritime traffic engine- 
ering was held in Swinoujście on the Baltic sea coast. 
It was organized by the Institute of Maritime Traffic Engi- 
neering, Maritime University of Szczecin. 

The vast program of the Conference comprised 66 pa- 
pers 6 of which were presented during the plenary ses- 
sion, and the remaining during 5 topical sessions. 


The following papers were presented 
during the plenary session : 


© Extreme parameters of ships intended for entering the 
port of Swinoujscie - by S. Gucma and W. Slaczka 
(Maritime University of Szczecin) 

> Mathematical model of ship motion in canals and locks 
by S. Zaikov, M. Lavrinovsky and V. Zaikov (State 
University for Water Communication, St. Petersburg) 

> Ship as intelligent machine - by R. Smierzchalski (Gdy- 
nia Maritime University) 

> Onnecessity of establishing a Polish Institute of Navi- 
gation following the example of similar institutions 
operating in neighbouring countries - by A. Weintrit 
(Gdynia Maritime University) 

> The competences and duties of the officers in the 
charge of a navigational watch in the face of coming 
into operation of the ship's integrated control systems 
by Z. Kopacz, W. Morgas and J. Urbanski (Naval Uni- 
versity, Gdynia) 

> A method for the improving of location accuracy of 
objects with the use of teledetection data - by J. Sanec- 
ki, A. Klewski, L. Cwojdzinski, K. Maj and P. Ka- 
mienski (Military Engineering Academy, Warsaw) 


The authors of the conference papers represented 15 
scientific research centres including one of Slovakia and 
one of Russia. The greatest contribution to the conference 
materials (36 papers) was made by the authors of Mariti- 
me University of Szczecin. The conference became a com- 
prehensive review of the interesting research projects con- 
cerning maritime traffic engineering, carried out in this 
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UNDERWATER TECHNOLOGY 


Kinematical control of motion 
of underwater vehicle 
in horizontal plane 


Jerzy Garus 
Naval University, Gdynia 


ABSTRACT 


In the paper presented is a method of designing a fuzzy-logic-based autopilot for control 
of horizontal motion of an unmanned underwater vehicle. The control system 5 synthesis 
was performed under the assumption that the vehicle can move with variable linear and 
angular velocities and the quantities possible to be measured are : position and orienta- 
tion of the vehicle in the inertial reference system. The task of the autopilot was to minimi- 
ze the mean squares of deviations from the motion trajectory given in the form of a broken 
line defined by the coordinates of successive turning points. To generate control signals 


three independent fuzzy PD controllers using the control principles based on the Mac Vicar-Whelan’s 
standard base, were applied. For the linguistic variables of each controller appropriate fuzzy sets were 
selected and linear membership functions of trapezoidal and triangular form were defined. The presented 
results of the simulation tests performed for the remotely operated underwater vehicle ,, Ukwiat”, with 
and without influence of disturbances resulting from sea current, confirm the proposed approach to be 


correct and effective. 


Key words : Underwater vehicle, autopilot, fuzzy logic 


INTRODUCTION 


An increasing interest has been given to underwater robo- 
tics in the last years. Currently, it is common to use unmanned 
underwater vehicles (UUVs) to accomplish such missions as : 
inspection of coastal and off-shore structures, cable mainte- 
nance, as well as hydrographical and biological surveys. In the 
military field they are employed in such tasks as surveillance, 
intelligence gathering, torpedo recovery and mine counter 
measures. The main benefits of usage of the UUVs can be the 
possibility of removing a man from the dangers of the under- 
sea environment, and of reduction in cost of exploration of 
underwater space. 

There are various categories of the UUVs. The most often 
used are remotely operated vehicles (ROVs). The ROV is usu- 
ally connected to a surface ship by a tether through which all 
communication is wired. The tether’s drag influences motion 
of the vehicle and may produce significant disturbances and 
energy loss. The ROV is equipped with a power transmission 
system and controlled only by thrusters. Simultaneously the 
spatial station-keeping or tracking of the underwater vehicle is 
a difficult task for a human operator, hence a supervisory con- 
trol has been developed to support its own intelligence and 
autonomy. 

Automatic control of underwater vehicles is a complex pro- 
blem due to their strongly coupled and highly nonlinear dyna- 
mic characteristics. Moreover, the dynamic characteristics can 
change depending on a chosen vehicle’s configuration suitable 


to its mission. In order to cope with the difficulties the control 
system should be flexible. An interesting review of classical 
and modern techniques adopted to control the dynamic beha- 
viour of unmanned underwater vehicles was presented in [1,3]. 
Nowadays fuzzy logic control systems have been successfully 
applied to a wide variety of mechanical systems [2,7,9]. The 
primary advantage of the fuzzy controllers is the possibility of 
easy incorporating heuristic knowledge of experts into a con- 
trol strategy. In this paper a fuzzy autopilot for tracking control 
of the ROV is described. 


TRACKING CONTROL 


The general motion of a marine vehicle of 6 degrees of free- 
dom (DOF) can be described by the following vectors [3] : 


n=lx,y,z,0,0,w]’ 
v= lu, v,w, p,q, r] (1) 
t =[X, Y, Z, K, M, N|" 

where : 


the position and orientation vector in the earth- 
-fixed frame 

the linear and angular velocity vector in the 
body-fixed frame 

the forces and moments acting on the vehicle in 
the body-fixed frame. 
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The nonlinear dynamic equations of motion 
can be expressed in matrix form as [3] : 


Mv + C(v)v + D(v)v +gm) =T 
n= J(n)v 


where : 


(2) 


M -inertia matrix (including added mass) 

C(v) -matrix of Coriolis and centripetal terms 
(including added mass) 

D(v) -hydrodynamic damping and lift matrix 

g(1)) -vector of gravitational forces and moments 

J(N) -velocity transformation matrix. 


Coordinate systems and tracking control 


For the conventional ROVs a basic motion is the move- 
ment in horizontal plane with some variation due to diving. 
They operate in crab-wise manner with 4 DOF and small roll 
and pitch angles which can be neglected during normal opera- 
tion. Therefore, it is purposeful to regard 3-dimensional mo- 
tion of the vehicle as the superposition of the motion in the 
horizontal plane and that in the vertical plane. Farther in the 
paper only the movement of the vehicle in the horizontal plane 
is considered. 

It is convenient to define three coordinate systems when 
analysing route tracking systems for the marine vehicle mo- 
ving in horizontal plane (Fig.1) [6] : 


* the global coordinate system OyXpYo (the earth-fixed fra- 
me) 

* the local coordinate system OXY (fixed to the body of the 
vehicle) 

* the reference coordinate system P;X;Y; (not fixed). 


xX. 


Oy Yi Yim Y Yo 


Fig. 1. Coordinate systems used to describe tracking control 


of the underwater vehicle moving in the horizontal plane: 
OpXo0Yo — global system, OXY — local system, P;X;Y; — reference system 


The main task of the designed tracking control system is to 
minimize the distance of attitude of the vehicle’s centre of gra- 
vity, d, to the desired trajectory under the following assump- 
tions : 


© the vehicle can move with variable linear velocities u, v, 
and angular velocity r 

> the vehicle’s position coordinates, x, y, and heading angle 
w are measurable 
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> the command signal T consists of three components : 
X and Y - forces along X- and Y-axis, respectively, 
and N - moment around Z-axis 
Ə a travel time is not given in advance, thus the navigation 
between two points is not constrained by time. 


The tracking autopilot has to provide both track-keeping 
and course-keeping capabilities. Hence, the autopilot should 
minimize the mean squares of deviations "d", from a desired 
track, and Aw deviations from a desired course: 


I= min > [d?(t)+ Ay? ()] 


where: 
d(t) = -sin(Ay)(x(t)- x; )+ 
+ cos(Ay)(y(t)—y;) 


(3) 


Yin Yi 


XiT Xi 


= 


@; =arctg 


Aw(t) 


wt) 
À 


- angle between the track reference 

line and vehicle’s centreline : y(t) - Q; 
- heading angle of the vehicle 
- constant coefficient. 


Each time the vehicle location x(t), y(t) 
at the instant t satisfies : 


ban = x(t)? t [Yia E yol < p’ 


where : 


(4) 


p - radius of circle of acceptance. 


The next waypoint should be selected on the basis of the 
reference coordinate system (e.g. Pi+1Xi+1 Yj+1) and the vehic- 
le’s position should be updated in compliance with the new 
reference coordinate system. 


Fuzzy control law 


A fuzzy proportional derivative controller, adopted from 
[2], working in the configuration presented in Fig.2, has been 
designed for the tracking control. 


Underwater 
vehicle 


Fuzzy controller 


Fig. 2. The assumed structure of the fuzzy controller 


The membership functions of fuzzy sets of the input varia- 
bles : the error signal e = naj - nj and the derived change of 
error Ae = nj - Nj-1 , aS well as the output variable : the com- 
mand signal 7; - for j € {1, 2, 6} - are shown in Fig.3. 


Values of unknown parameters: Xe, XAe, Xs and xy used in 
computer simulations are given in Tab.1. Evaluation of the 
parameters can be done by means of many optimisation tech- 
niques, classical or modern ones, e.g. Genetic Algorithms [4,5]. 


Tab. 1. The assumed membership function parameters 


Controller 


Zen 


Degree of membership 
oO j=) © t 
N > ON oo 


Degree of membership 
> > 2 2> 
> b BR ^ o 


-1 -Xhe 0 Xis 1 
Derived change of error Ae 


Degree of membership 


Command signals X,Y,NV 


Fig. 3. The assumed membership functions for fuzzy sets of : 
the error e, derived change of error Ae, and command signals X, YN. 
Notation : N - negative, Z - zero, P - positive, 
S - small, M - medium and B - big 
The chosen control rules, taken from the Mac Vicar- 
-Whelan’s standard base of rules [8], are given in Tab.2. 


Tab. 2. The assumed base of rules 


Error signal e 


Derived 
change of 
error Ae 


Command signals X, Y,N 


SIMULATIONS 


A simulation study of tracking control has been performed 
for an underwater vehicle ,,Ukwiat’” designed and built by 
Gdansk University of Technology for the Polish Navy. The 
ROV is an open, 1.5 m long frame robot controllable in 4 DOF, 
fitted with a propulsion system consisted of 6 thrusters. 

The vehicle can move in horizontal plane by using four thru- 
sters. Every thruster can generate thrust force up to +250 N. 
It assures its speed up to u = +1.2 m/s and v = +0.6 m/s in X and 
Y direction, respectively. The autopilot in question consists of 
3 independent controllers producing the command signals X, Y 
and N calculated on the basis of the proposed fuzzy law (under 
the constraints : |X| < 500 N, |Y|< 150 N and |N|<50 Nm). 


The structure of the proposed control system 
is presented in Fig.4. 


Numerical simulations have been made to confirm validi- 
ty of the proposed control algorithm under the following as- 
sumptions: 


* the vehicle has to follow the desired route beginning from 
the position and orientation point : (x = 10 m, y = 10 m, 
yw = 0°), passing the target waypoints : (10 m, 90 m, 90°), 
(30 m, 90 m, 0°), (30 m, 10 m, 270°), (60 m, 10 m, 0°) and 
ending at the point : (60 m, 90 m, 90°) 

æ the turning point is reached when the vehicle operates wi- 
thin the 1,5 m radius of circle of acceptance, p 

* the initial conditions are the same. 


Flight Planner 


Autopilot 


trajectory | trajectory | course 
in X-axis | in Y-axis 


Underwater 
Vehicle 


Fig. 4. The block diagram of the track-keeping system 
Notation : index d - desired value 


The tracking control simulation results and the courses of 
command signals for no added environmental disturbances are 
shown in Fig.5. The real route and position of the vehicle al- 
most coincides with those desired. Also the quality of course- 
-keeping control is satisfactory (course deviations close to zero). 
In Fig.6 illustrated is an influence of sea current disturbances 
on the vehicle’s route and course. A clear difference between 
the desired and real tracking curves is there observed. Altho- 
ugh the errors of position and course are much bigger than in 
the previous case the autopilot is able to cope with the external 
disturbances and to reach the turning points with the comman- 
ded orientation. It should be noted that the last simulations were 
performed for the sea current speed increased to 0.35 m/s. 
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Fig. 5. The vehicle 5 tracking curve and course deviations from the desired position and course, 
for the command signals without interaction of environmental disturbance 
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Fig. 6. The vehicles tracking curve and course deviations from the desired position and course, 


for the command signals under influence of the sea current of 0.35 m/s speed, and 135° direction angle 


desired 


Notation : 
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real 


This is comparable with the average speed of the vehicle. As it 
was expected, the guidance was very sensitive to the ratio of the 
vehicle speed and the current speed. It was observed that this 
ratio should not be smaller than 2, which is also confirmed in [6]. 


CONCLUSIONS 


“+ In this paper the waypoint-tracking autopilot using fuzzy 
control, intended for underwater vehicles, has been descri- 
bed. The nonlinear model of the ROV ,,Ukwiat’” was ap- 
plied to carry out computer simulations. 

«* The simulation results obtained by using the control sys- 
tem design method based on three decoupling fuzzy con- 
trollers, showed the presented autopilot to be simple and 
useful for practical use. 

* The main advantage of the proposed solution is its flexibi- 
lity with regard to the vehicle’s dynamic model, and its high 
performance at relatively large disturbances resulting from 
sea current. 


NOMENCLATURE 

C(v) - matrix of Coriolis and centripetal terms (including added 
mass) 

d - distance of attitude of the vehicle's centre of gravity 


from the desired trajectory 


DOF _ - degress of freedom 

D(v) - hydrodynamic damping and lift matrix 

e - error signal 

gm) -vector of gravitational forces and moments 
J - quality index 

J(n) -velocity transformation matrix 

K - moment around X- axis 

M - inertia matrix (including added mass) 

M - moment around Y- axis 

N - moment around Z- axis 

OXY -local coordinate system (fixed to the body of the vehicle) 


OpXoYo - global coordinate system (the earth-fixed frame, inertial) 
p - angular velocity around X- axis 
P;X;Y; -reference coordinate system (not fixed, connected 
with trajectory of motion) 
- angular velocity around Y- axis 
- angular velocity around Z- axis 
- time 
- linear velocity along X-axis 
- linear velocity along Y- axis 
- linear velocity along Z- axis 
- vehicle's position coordinate along Xo - axis 
es Xe - coordinates of membership functions of e and Ae, 
respectively, (along axis of abscissae) 
Xs, Xm -Coordinates of membership functions of T;, (along axis of 
abscissae) 
X - force along X- axis 
y - vehicle's position coordinate along Yo - axis 
Y - force along Y- axis 
Z - vehicle's position coordinate along Zo - axis 
Z - force along Z- axis 
A 
A 


g<¢e77 0 


no o 


e - derived change of error 
y - change of heading angle (deviation from a desired 
course) 
n - vehicle's position and orientation vector in the earth- 
-fixed frame 
0 - vehicle's pitch angle 
À - constant coefficient 
v - linear and angular velocity vector in the body-fixed 


frame 

p - radius of circle of acceptance 

T - forces and moments acting on the vehicle in the body- 
-fixed frame 

Q; - direction angle of i-th segment of trajectory 

o - vehicle's roll angle 

y - vehicle's yaw (heading) angle 
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FOREIGN 


A Hundred-Year Jubilee 


On 7 July of the year 1903 — in accordance with the 
record in the commemorative book — the research establi- 
shment for water engineering and shipbuilding VWS (Ver- 
suchanstalt fuer Wasserbau und Schiffbau) officially com- 
menced its activity in Berlin. 


In the course of time it could boast of outstanding - in 
the worldwide scale - scientific and technical achievements 
in the domain of inland waterways, harbour engineering 
and shipbuilding. 


On this special occasion the German society of naval 
architects and marine engineers STG (Schiffbautechnische 
Gesellschaft) held its general assembly on 19 - 22 Novem- 
ber 2003 just in Berlin. 


Its agenda contained , apart from its organizational part 
dealing with STG matters, the special session devoted to 
the achievements of the VWS, as well as the scientific con- 
ference carried out within 5 topical sessions. 


Scientific workers of the Faculty of Ocean Engineering 
and Ship Technology, Gdansk University of Technology 
have had multi-topical contacts both with the VWS and 
STG. Their representative, D.Sc.Eng. Edmund Brzoska 


C part in these events. D 
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Faculty of Ocean Engineering 
and Ship Technology 


Gdansk University of Technology 


New Missions And Visions 
Towards European 
Research Area 


The Faculty of Ocean Engineering and Ship Technology 


of Technical University of Gdansk has recently revised the 
objective of activities according to European Research Area 
assumptions. In this regard the new additional missions of 
Faculty are set: 


> 


> 


32 


to enhance the productivity and effectiveness of the ma- 
ritime sector of Poland as a member of EU 

to promote fast and efficient integration of Polish mari- 
ne industry into EU and ERA (European Research Area) 
to provide the efficient structure for exchange of rele- 
vant knowledge and information within EU through or- 
ganisation of international conferences, workshops and 
schools 

to promote a co-ordinated research effort in the mariti- 
me sector of the European Research Area, in particular 
to assist in creation of consortia participating in the Eu- 
ropean Framework Programmes 

to organise and support the short- and long-term exchan- 
ge visits of researchers from and to EU 

to attract the promising young people into the research 
and production enterprises of the European maritime tech- 
nology sector 


> 


> 


to organise and execute continuous training of highly qua- 
lified researchers through D.Sc. studies 

to organise and execute continuous training of designers 
in the most advanced and sophisticated design techniqu- 
es through post-graduate courses 

to organise and execute continuous training of shipyard 
and co-operating industry staff in implementation of the 
advanced production techniques 

to organise and execute continuous training of the ship- 
-owners and relevant state administration personnel in 
the problems of safe and reliable operation of the sea 
transport and offshore industry 

to provide expertise and highly qualified professional 
consultations to small and medium enterprises of the ma- 
ritime sector of NAS (Newly Associated States). 

to create the efficient and well-organised structure for 
systematic collection, generation, exchange and disse- 
mination of advanced and highly specialised knowledge 
between research organisations, educational institutions, 
shipyards and associated industries, ship-owners and state 
administration 

to initiate and continuously support a co-ordinated effort 
towards design, production and operation of durable, re- 
liable, economic and technologically advanced ships. 
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